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Abstract 

If non-vanishing chemical potentials are assigned to chiral fermions, then a Chern-Simons 
term is induced for the corresponding gauge fields. In thermal equilibrium anomalous pro- 
cesses adjust the chemical potentials such that the coefficient of the Chern-Simons term 
vanishes, but it has been argued that there are non-equilibrium epochs in cosmology where 
this is not the case and that, consequently, certain fermionic number densities and large-scale 
(hypermagnetic) field strengths get coupled to each other. We generalise the Chern-Simons 
term to a real-time situation relevant for dynamical considerations, by deriving the anomalous 
Hard Thermal Loop effective action for the hypermagnetic fields, write down the correspond- 
ing equations of motion, and discuss some exponentially growing solutions thereof. 
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1. Introduction 



It was realised long ago that if chiral fermions are assigned a non-zero chemical potential, 
then a Chern-Simons term appears to be induced for the corresponding gauge fields T!. This 
statement is not without ambiguities, however [2]. In fact a non-Abelian Chern-Simons term 
transforms non-trivially under large gauge transformations, and therefore the induced theory 
is well-defined only for certain imaginary values of the chemical potential Formally, the 
reason for these problems is that the chiral charge is not conserved because of the axial 
anomaly 1|, so that strictly speaking no chemical potential should be assigned to it. 

In the standard electroweak theory, chiral fermions couple not only to non-Abelian gauge 
fields, but also to the Abelian hypercharge fields. Hence, at the high temperatures where 
the electroweak symmetry is restored, a Chern-Simons term appears to be induced for them 
as well. An Abelian Chern-Simons action does not have the same topological properties as 
the non-Abelian one, but on the other hand it is gauge-invariant, and could thus conceivably 
have a more direct physical significance than its non-Abelian counterpart. 

Indeed, there have been a number of suggestions for possible roles that the Abelian hy- 
permagnetic Chern-Simons term might play in cosmology. One of them is related to the 
observation that right-handed electrons, which do not take part in weak interactions and 
also have a very small Yukawa coupling, are practically decoupled from the thermal ensemble 
above temperatures of about 10 TeV 131. If they come with a non- vanishing net density, which 
can be described by a chemical potential, then a hypermagnetic Chern-Simons term gets in- 
duced. It has been argued that this leads to an instability and to the subsequent generation 
of large-scale hypermagnetic fields jH]. While it is believed that any length scales related to 
physics within the horizon of this epoch are too small to act as seeds for the currently observed 
galactic magnetic fields |3|, such fields could have other physical consequences, for instance 
affecting the properties of the electroweak phase transition p|, the sphaleron energy [Hj, and 
electroweak baryogenesis ^01 (in suitable extensions of the Standard Model). 

Another possible role acts in the opposite direction. Suppose that there exist primordial 
hypermagnetic fields as a result for instance of some inflationary dynamics. Then anomalous 
processes could convert some of these fields to lepton and eventually to baryon number, 
resulting possibly in the existence of matter and antimatter domains which could lead 
to nucleosynthesis taking place in a corresponding environment jl2l . 

Whichever of these physics effects is realised, a central ingredient is always the presence of 
a hypermagnetic Chern-Simons term induced by a fermionic chemical potential. As argued 
in Refs. the Chern-Simons term leads to an additional "anomalous" term in the 
magnetohydro dynamic equations that govern the evolution of the hypermagnetic fields. The 
goal of this paper is to attempt a field theoretic derivation of the equations of motion for 
the hypermagnetic fields in the presence of fermionic chemical potentials. The appropriate 
framework is that of the Hard Thermal Loop (HTL) effective theories |13lll4j . By integrating 
out the "hard" fermions, with energies of the order of the temperature T, we derive the 
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effective action for the "soft" gauge fields, with wave vectors p and Minkowskian frequencies 
uj much smaher than the temperature, |p[, \uj\ <ti 2ttT. The standard static Chern-Simons 
term used for instance in the considerations of Refs. OCI] is recovered if we make the further 
approximation \lo\ <C |p| (which indeed appears to be well justified in practice). 

The plan of this note is the following. In Sec. [21 we recapitulate the static gauge field 
effective action at high temperatures, and in Sec. generalise it to the non-static situation. 
We solve the resulting equations of motion for a simple case in Sec. 01 and conclude in Sec. El 



2. Anomalous effective action in the static limit 

Let us consider the standard electroweak theory at temperatures T above a few hundred 
GeV. Let be the SU(2)l and U(1)y gauge fields, C^^^Ff^i, the corresponding field 

strength tensors, and g, g' the gauge couplings. The covariant derivative reads = — 
igT'^A'^ + ig'^B^ , where are Hermitean generators normalised as Tr [T'^T''] = 5°"^ /2, and 
Y is the hypercharge quantum number. With these conventions, the gauge field part of the 
dimensionally reduced Euclidean Lagrangian |E1 takes the form 



fE+ ^ GIliyG^lu + ^ F,,uF, 



,,u- ,,u + \ ni\A%A% + ^ m'iBoBo + ijEBo + 



+ c^ncs + c'^in'cs + ... 



(2.1) 



where the "anomalous" Chern-Simons densities read 
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and infinitely many higher dimensional operators have been suppressed. The term is a 
(field-independent) "unit operator". The relative parametric accuracy that can be reached 
with dimensionally reduced effective theories of this kind has been analysed in Ref. HW. 

There are a number of matching coefficients appearing in Eq. ()2.1() . Denoting by no the 
number of generations, by the number of fundamental scalar doublets, by the common 
chemical potential of the quarks, and by /U^. (/i/?,) the chemical potential of the i^^ left-handed 
(right-handed) lepton generation, the leading-order expressions read 
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Some higher-order corrections can be found in Ref. |17j . 

Now, because of the axial anomaly, the rate of baryon plus lepton number violation, 
din \B + L\/dt a -(13nG/4)(25.4± 2.0)a^r HH-Hl],-^ is significantly larger than the expan- 
sion rate of the Universe, ~ T^/mpi, for 10^ GeV <r< 10^^ GeV, so that the anomalous 
processes are perfectly in thermal equilibrium ::23_. Therefore, the corresponding chemical 
potential should be set to zero: 



(2.10) 



i=l 



In other words, the coefficient ce in Eq. (|2.8I) vanishes. 

On the other hand, the coefficient in Eq. (|2.9() does not vanish, provided that nn- ^ ^l.. 
Such a situation can arise if chirality flipping processes, mediated by the Yukawa couplings, 
are out of equilibrium |^ , and we will assume this to be the case in the following. Formally, 
this can be reached by setting the electron Yukawa coupling to zero. 

The naive conversion of Eq. ^2.1\ to Minkowski spacetime goes simply through the analytic 
continuation 
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and the Minkowskian action is then given by Sm = /dt d^x jCj\/. The resulting theory is 
gauge invariant only in static gauge transformations, however, and thus cannot be the full 
truth. In the next Section, we recall how a more precise theory can be obtained. 



3. Anomalous Hard Thermal Loop effective action 

Suppose that we stay for a further moment in the static limit, and consider what kind of 
higher order operators could appear in Eq. 1)2.1(1 . From the point of view of the original 
four-dimensional theory, some of these operators arise from a gradient expansion in spatial 
derivatives. Given that the scale that has been integrated out to obtain Eq. (|2.1() is the 

^The number 25.4 is in fact the value of a function containing terms like ln(l/Quj), at the physical oiw 
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"hard" scale ~ InT, these operators are necessarily suppressed by (at least) 0(| Vp/(27rT)^) 
with respect to the ones that have been kept in Eq. ()2.1() . 

Now, one might expect that the same is true for temporal derivatives: maybe their effects 
are also suppressed by / {2ttT)'^7 This is not the case! As is well-known from Hard Thermal 
Loop considerations j^ElEI, time-dependence is suppressed with respect to the static limit 
only by Od^ol/lVl), and is in general of order unity. 

Let us proceed with the explicit computation. We employ the Matsubara formalism, fol- 
lowed by analytic continuation. The chemical potential corresponds in momentum space to 
shifting fermionic Matsubara frequencies w„ as ojn w„ + i/i. We use the (sometimes im- 
plicit) notation that summing over Lorentz indices which are both down implies the use of 
Euclidean metric. Capital momenta (P, Q, R) are assumed Euclidean. 

With these conventions, for any given left-handed fermion with hypercharge Y, the anoma- 
lous part of the hypermagnetic Euclidean action in momentum space reads 



6Se 



Ig'^Y'^ 6q+rB^{R)B,{Q)T^,{Q) , (3.1) 
4 Tq,r 



r,AQ) - (P + Q)2p2 ' (3-2) 

^^u{Po,p) ^ ^Tr[7^(f +^)7.f75] , (3.3) 

where f] denotes the usual Matsubara sum-integral (bosonic for Q,R and fermionic for P). 
For right-handed fermions the overall sign is opposite. The Euclidean 7-matrices here are 
Hermitean, with {7^,71/} = 25^^, and we have defined 75 = 70717273- It is perhaps appro- 
priate to mention that it is possible to have an "anomalous" term with two gauge field legs 
only, since the chemical potential effectively acts as a third leg. 

To evaluate the sum over uJn in Eq. (|3.2|) . we note that there are only single poles in po 
and that we can thus use the contour formula 

n odd Im z < Im z > 

where the sums are over the poles z of f{z). The first term is independent of T and /i; we 
ignore this zero-temperature vacuum part here. The latter two terms are ultraviolet and 
infrared finite, and require no regularization. Thus we can work in exactly four dimensions 
(as already implicitly assumed above), whereby Tr [757/x7i^7a7/3] = 4e^i,Q/3, with £0123 = +1- 
Picking up the four poles; shifting integration variables in three of them^ as p ^ — p, 
p ^ p - q, p ^ -p - q; and employing the facts that A^^{-po,-p) = -A^^{po,p), 
^fiuiPo - Qo, P - q) = ^fiuiPo, P), ^fiu{-Po - go, -P - q) = -^fj.u{po, P); the thermal part 

^Shifting integration variables is safe provided that each individual term is finite, as is the case here. We 
have however checked the outcome also by not carrying out any shifts but just expanding in Q/cup. 
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of Eq. (|3.2|1 becomes 



Q + 



Q 



(3.5) 



where we have denoted Jp = /d^p/(27r)^, N-{ujp) = nY{u}p — fi) — nF(a;p + ^), nF(wp) = 
l/(e^'^p _j_ -^^^ yE ^ (^—i^p^/ujp), ujp = |p|, and ■ Q = vj^Q^. Note that no approximations 
have been made so far for the thermal part. 

The next step is to carry out a smah couphng expansion. In other words, we look for 
the leading term in the expansion in small \Q\/ujp, where parametrically (after the analytic 
continuation to follow presently) Q is a soft scale, \Q\ <max{gT,gfi), while the integration 
variable gets its contributions from the hard scales, ujp ~ max(T, fj,). The leading term in the 
expansion of the denominators in Eq. H3.5|) obviously cancels, but the next-to-leading term 



is non-vanishing, and multiplied by A^iy(— iup, p) 



^p ^fia, 



upQaVp ■ Thus, we obtain 
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At this point we carry out analytic continuation to Minkowski spacetime. Any Euclidean 
Lorentz-vector can be written as = Aq,/?/*^^, with K^p = diag(— i, —1, —1, —1). Further- 
more, as mentioned in Eq. 1)2. there is an overall minus-sign between Ce and Cm- We 
also introduce an angular integration = /dr2^,/47r, where = (l,t'*), v^v^ = 0, and the 
integral is over the directions O^, of u*, and note that the radial integration in Eq. ()3.6j) can 
be carried out exactly, with the result 
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Summing over all fermions, adding the known non-anomalous bosonic terms ^1], and going 
furthermore to j;-space, we arrive at our final action: 
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.Ff^" and = d-d 



d'^d,, - The coefficients m'^,j'^, c'^ are the same as in 



Eqs. (Eini), (EIZI, and (EIHI)-^ The integrals J^{...) are collected in Appendix A. 

It is the last term in Eq. ()3.9() which is the new one. For more symmetry, one could replace 

since the extra term —[d^/{v - d)\v ■ B thus introduced vanishes 
after partial integration and neglectance of surface terms. 

^We note that a single chiral fermion contributes to the parameters as j'^ = g'fiY{T^ + /i^/7r^)/12, m'^ = 
g'^Y^{T^/3 + tJ-'^/Tv'^)/8, c'e = -fj.Y^H/2, where H ^ +1 (-1) for right-handed (left-handed) fermions. 
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Note that Eq. ()3.9|) is gauge invariant in time-dependent gauge transformations, unlike 
Eq. ()2.1() . Nevertheless, by employing the identities in Eqs. HA.11|) . (|A.12|) . it is straightfor- 
ward to rewrite Sm in a form where it is obvious that the hypercharge part of Eq. (|2.1|) is 
recovered in the static hmit, up to corrections of order 0(|5o|/|V|). 

For completeness, we remark that although we have shown in Eq. (|3.9|) only the terms with 
the three largest coefficients, j'^ ~ 0{^T'^),ni''^ ~ 0{T'^), and c'^ ~ C>{n), some additional 
operators are known as well: there is a non-anomalous but charge-conjugation violating 
"cubic" operator with a coefficient 0{n) and it appears that even "quartic" operators 
with a coefficient 0(1) could be added, by reformulating the HTL action in terms of a certain 
matrix- valued kinetic theory j2Hj . 

Since the action in Eqs. 1)3. 8() . H3.9() is quadratic in the fields, its contents can equivalently 
be expressed through equations of motion. It is useful to express the equations of motion in 

XT2I). We obtain 



a form following from the identities in Eqs. (jA.ll 

/2 VuV'^ 
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The integrals have the properties v°' = (5°o and (1 — 6°'o) v°' /v • 9 oc 5" (cf. Appendix A), 
which imply that the right-hand side is divergenceless (or transverse) with respect to d"^. 

We should like to stress already at this point that there are circumstances where higher 
order corrections to Eq. ()3.10|1 can become important. Let us recall, to start with, that 
standard HTL structures like the second term in Eq. 1)3. 1U() can be reproduced by classical 
kinetic theory, or "Vlasov equations" , of the type 



p-da^-^p''F^fsd^]f{x,p) 
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by solving them perturbatively to second order in g' with the "initial condition" that for 
g' —>■ the solution reads f{x,p) = f{o){x,p) = A'K5{p^)n-p{\pQ\ — sign(po)M)- One may then 
expect higher order interactions to generate a Boltzmann type collision term on the right- 
hand side of Eq. 1)3.11(1 . Indeed, effective theories which such a structure have recently been 
analysed in Ref. and the collision term does turn out to be important in many contexts. 
We will return to this issue presently. 

Finally, the question arises whether the last term in Eq. ()3.10() can also be expressed 
through Vlasov equations. In principle this indeed is the case: for instance, we can "by brute 
force" generalise f{x,p) to a two-index Lorentz tensor, and the Vlasov equations to 



P ■ dx 
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(3.14) 
(3.15) 



6 



where H = +1 (—1) for right-handed (left-handed) fermions. Solving this set perturbatively 
to second order in g' , it is easy to check that Eq. (|3.1U|) is reproduced (with a single chiral 
fermion contributing as specified in footnote 3). Nevertheless the new parts in these equations 
are not particularly satisfactory: Eq. (|3.13|) generates a distribution which contains an ugly 
S'ijP') and is not gauge invariant. While both problems disappear after the construction of 
the current in Eq. (|3.14|) . it should be possible to find a more compelling formulation which 
also has a physically understandable interpretation. 

4. On the growth rate of instabihties 

In order to illustrate the effects that may originate from the non-local HTL structures, we 
inspect the spatial part of Eq. (|3.1fl|) (i.e. i/ = 1, 2, 3). We employ a class of gauges where Bq 
is constant, so that F^q = —doBk, where k is a spatial index. 

We start by considering the static limit. Then Eq. (|3.1())) becomes {k = 1, 2, 3) 

[V^4i - dkdi]Bi{x) = CEY^ekijFijix) . (4.1) 

Going into momentum space [Bi^{x) = i?fc(q) exp(— zq • x)], and choosing a frame where 
q = (0,0, q^), we obtain equations for the transverse components Bi,B2. They have a non- 
trivial solution provided that 

lql = ±|fcs|, 9cs = ^. (4.2) 

Going back to configuration space, the solution reads 

Bi{x^) = C cos[gcs(x' - xl)] , 52(x=^) = C Mqcsix^ - xg)] , (4.3) 

where C,Xq are constants. This is nothing but the so-called (static) Chern-Simons wave j29j . 

Consider then the dynamical situation. We transform Eq. (|3.1U() to Fourier space with 
respect to space coordinates but keep the time coordinate in configuration space. Building 
on Eq. 1)4. we look for transverse modes which satisfy the eigenvalue equations {k' = 1,2) 

c^Bk' + iqcsek'ijQiBj = XBk' ■ (4.4) 

Two non-trivial solutions exist, with the eigenvalues 

A = q2±fcslql. (4.5) 

Thus, for long wavelengths, |q| < \qcs\j there exists an unstable branch with A < 0. 

How fast do the modes with A < grow? To find out, we assume that the time evolution 
is very slow, \df \ <C |A|, and justify this assumption a posteriori. In this situation, we can 
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approximate the complicated HTL structures through the leading terms in time derivatives, 
and the problem becomes tractable. 

Inspecting the integrals in Appendix A, it can be seen that both the 2nd and the 3rd 
term on the right-hand side of Eq. (|.S.in|) have a term linear in time derivatives. For a 
solution of Eq. (|4.4|) . however, the third term gives a contribution suppressed by C[(q^ — 
A)/m^] ~ 0[(7^g/m^] with respect to the second term. This is very small for weak coupling, 
~ I max(r^, //^)7r^, and can safely be ignored. Thus the effect comes from the 2nd term, 
and we obtain 

- \Bk' ~ -^^IdtBy . (4.6) 
To summarise, modes with A < grow exponentially. By ~ exp(rt), with the rate 

P _ 4|Aq| ql^ 

r~ — 72 ~ Z;72lql • (4.7) 



E 



Given that q^^/m''^ <C 1, the assumption of slow growth (F^ ^ q^) is indeed justified. 

It is important to realise at this point, though, that higher order corrections may give 
large contributions on the right-hand side of Eq. (|4.6j) . as already mentioned above. In fact, 
it could happen that summing an infinite set of higher order loop contributions effectively 
shields the scale |q| in the denominator of 7rmg/4|q| by a constant ~ g"^\n{\/ g')T , whereby 
the right-hand side goes over into a'dtB^i, where a' is the hyperelectric conductivity, a' ~ 
T/g''^hi{l/g') 30 . Another way to think of the issue"^ is that interactions may generate 
a "thermal width" Fth of order Fth ~ g"^\n{l/ g')T for the hard on-shell particles and 
that, if qo ^ qo + «Fth and |q| ^ Fth in Eq. (|A.10|) . then the static limit would become 
— i/Fth rather than — i7r/2|q|, amounting to the same shielding. Since |q| ~ qcs ~ g''^^J' and 
Tth ^ ln(l/(7')r, we should formally assume |q| ^ Fth so that the shielding is irrelevant, 
but since typically ^ <C g'^T in cosmology, the formal hierarchy may get reversed so that we 
indeed find ourselves in the situation |q| ^ Fth- In fact, inserting numerical values relevant 
for the Standard Model, it appears that having | Y^'^^iit^ii ~ f^njl/'^ ^0.5 already brings us 
to the reversed situation. 

The estimates for the growth rate of hypermagnetic fields that were presented in Refs. |H|lllj 
were based on the contribution of the hyperelectric conductivity a', rather than Eq. 1)4. 6|) . and 
should thus be correct under the phenomenologically relevant circumstances fj, <^ T. (The 
expansion of the Universe as well as the time-dependence of the chemical potentials, due to 
the backreaction via the hypermagnetic part of the anomaly equation as well as the chirality 
flipping processes induced by the electron Yukawa coupling, were also taken into account.) 
We note that the corresponding growth rate is much larger than Eq. 1)4. 7|) . since Fth ^ |q|- 

*We stress that the discussion here is only qualitative in nature, and omits important points. 
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5. Conclusions 



We have addressed here the couphng between hypermagnetic fields and fermionic chemical 
potentials in the standard electroweak theory at high temperatures. This problem has phe- 
nomenological relevance in cosmology, provided that a lepton asymmetry and/or primordial 
hypermagnetic fields exist at temperatures above the electroweak phase transition, but is 
also related to some intriguing theoretical issues, such as that the coupling discussed seems 
to allow for a sharp distinction between the high-temperature and low-temperature phases 
of the electroweak theory |32j . 

Concretely, we have generalised the standard Abelian Chern-Simons term to an apparently 
Lorentz invariant form, which can be added to the Hard Thermal Loop action describing 
the real-time dynamics of the hypermagnetic fields (Eq. ()3.9() 1. We have also analysed the 
unstable exponentially growing solutions that the resulting equations of motion have. Our 
conclusion is that for such solutions, the deviation of the anomalous term from its standard 
static form is in fact insignificant in practice (cf. the discussion preceding Eq. (|4.6|) ). so that 
the ignoring of this deviation in previous studies appears well justified in retrospective. 

The problem with the Hard Thermal Loop equations of motion is that higher order cor- 
rections to their non-anomalous part turn out to be very important for the small values of 
chemical potentials that are assumed to appear in cosmology, fi <^ T. (The wave vectors 
of the growing modes are proportional to differences of chemical potentials and thus very 
much smaller than the temperature in this situation.) In particular, the fact that a finite 
conductivity is expected to be generated through summing infinitely many high order loop 
corrections modifies the growth rate of the unstable solutions significantly. (The existence of 
unstable modes is not affected.) Nevertheless, it seems to us that the conclusion mentioned 
above, namely that it is safe to use the static limit of the Chern-Simons term under phe- 
nomenologically relevant circumstances, continues to be valid. To be sure, it would of course 
be interesting to develop a numerical framework where both the Hard Thermal Loop effects 
discussed in this paper, and higher order corrections such as conductivity, can be incorporated 
simultaneously, to allow for a more precise study of growing hypermagnetic fields. 

Finally, we remark that the growing hypermagnetic fields have a form which generates a 
non-zero value for the hypermagnetic topological charge density, ~ g''^F^uFfMu- This means 
that fermion number densities and, consequently, chemical potentials, evolve as dictated by 
the anomaly equation, in a way which stops the growth at some point. These processes have 
been analysed in Refs. |6l ITT]. 
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Appendix A. Basic integrals 

Although well-known |24l 1131 114j , we collect here the basic velocity integrals needed in this 
paper. Assuming implicitly that the frequency is replaced everywhere through qq — > qQ + iO^ , 
as is relevant for retarded Green's functions, the integrals read = 1,2,3) 
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where = (l,w*), q = {q^,c[), our metric convention is (H ), and 
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(A.4) 

(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 



(A.IO) 



Integrals with higher powers of f • g in the denominator can be obtained through the partial 
derivatives d/dq^. The following identities (which can be derived by certain partial integra- 
tions, or by explicit inspection) are often very useful: 



Jv (v ■ qY JvV-q 



V (v-q) 



{v ■ qY 



V V ■ q 



(A.ll) 
(A.12) 



Here 



q^ila, and /, J are arbitrary Lorentz vectors. 
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